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In the present work, we calculate the conductance through a single quantum dot weakly coupled 
to metallic contacts. We use the spin-1/2 Anderson model to describe the quantum dot, while 
considering a finite Coulomb repulsion. We solve the interacting system using the non-crossing- 
approximation (NCA) and the one-crossing approximation (OCA). We obtain the linear response 
conductance as a function of temperature and energy position of the localized level. From the com- 
parison of both approximations we extract the role of the vertex corrections, which are introduced 
in the OCA calculations and neglected in the NCA scheme. As a function of the energy position, 
we observe that the diagrams omitted within NCA are really important for appropriately describing 
transport phenomena in Kondo systems as well as in the mixed valence regime. On the other hand, 
as a function of temperature, the corrections introduced by OCA partly recover the universal scaling 
properties known from numerical approaches such as the Numerical Renormalization Group(NRG). 



I. INTRODUCTION 

Large- N expansions, with representing the angular 
momentum degeneracy, are commonly used for solving 
the Anderson Impurity model (AIM) and to study the 
Kondo physics. In particular, the so-called non-crossing 
approximation (NCA) in its infinite Coulomb repulsion 
limit captures the formation of the Kondo resonance at 
finite temperatures^ When compared with the numeri- 
cal renormalization group (NRG), the NCA scheme also 
provides the correct Kondo temperature {Tk)~ This suc- 
cessful match results from the fact that NCA collects, in a 
self-consistent way, all diagrams up to order 1/N. How- 
ever, it is well known that the Fermi liquid properties 
are not properly described. To deal with this deficiency, 
infinite resummations of spin flip terms are required 

When a finite value of Coulomb repulsion U is allowed, 
the NCA gives a severely underestimated value of Tk& 
The leading crossing diagrams that restore the proper 
energy scale were introduced in the early work by Pr- 
uschke and Grewe within the framework of the enhanced- 
NCA (ENCA) also called One Crossing Approximation 
(OCA)ii Unlike NRG and Quantum Monte Carlo, which 
are actually more robust techniques, the OCA method 
can be extended to multi orbital systems without much 
numerical effort. It also works on the real axis and can go 
to temperatures far below the Kondo one. These features 
make OCA a powerful impurity solver when combined 
with electronic structure calculations within the context 
of the Dynamical Mean Field Theory (DMFT)4 

Due to the interesting advantages mentioned above, it 
is important to check the role of vertex corrections on 
diagrammatic techniques in the calculations of different 
properties. Recently, the influence of vertex functions in- 
corporated by OCA has already been studied for general 
lattice problems^ multi-orbital Anderson models^ and 
dynamic susceptibilities of the AIMjii among others. 

In this work, we study the effect of vertex corrections 



on transport properties, in particular transport through 
a quantum dot (QD) weakly coupled to metallic con- 
tacts. The NCA results for the conductance, in the 
linear-response regime, were previously analyzed by Ger- 
ace et al.*£ The aim of our contribution is to compare 
the equilibrium conductance, as a function of tempera- 
ture and gate voltage, as calculated by NCA and OCA 
methods. Furthermore, we analyze the scaling properties 
of the conductance by comparison with the empirical for- 
mula extracted from NRG calculations. 

We observe that the diagrams omitted within NCA 
are really important for appropriately describing trans- 
port phenomena. Furthermore, the vertex corrections in- 
troduced by OCA partially recover the universal scaling 
behavior of the conductance as a function of temperature. 

The paper is organized as follows. In Sec. [II] we intro- 
duce the model and the calculation methods. In Sec. IIIII 
we present and discuss numerical calculations. Finally, 
in Sec. HVI some conclusions are drawn. 

II. MODEL AND METHOD 

To describe the system in which the QD is coupled to 
two leads we use the spin —1/2 Anderson Hamiltonian 

H = ^ £kv a J kl , a c kva + ^ E <mda + U n d1 -n dir + 

kva (J 
kua 

where the operator c kv „ represents a conduction electron 
with momentum k and spin a in lead v, where v = L, R 
labels the left and right leads; the operator d\ stands for 
an electron in the dot and nda is the number operator for 
a given spin in the QD. The parameters Ed and U rep- 
resent the energy of a single electron and the Coulomb 
repulsion in the dot, respectively. The coupling param- 
eter and hybridization functions between leads and QD 
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are given by V kv and A v (u) = ir J2 Ka V* u S(oj - e kva ), 
respectively. 

Our starting point is an auxiliary particle representa- 
tion of the Hamiltonian given byeq. In this repre- 
sentation, the physical operator is given by the following 
combination of the auxiliary particles, d\ = f^b+aa^ f- a . 
Operators b, f a and a label vacuum, single, and doubly 
occupied states, respectively. Within the OCA scheme, 
the self-energies of the auxiliary particles and the vertex 
functions are obtained self-consistently. The set of OCA 
equations are sketched in Fig. [TJ 7 



*i ~ w\aa# + tera# 

vv \ = 




FIG. 1. Sketch of the self-consistency OCA scheme for self- 
energies and vertex corrections of the auxiliary particles. The 
full, wiggly, dashed, and curly lines stand for pseudofermion, 
empty boson, conduction electron, and doubly occupied boson 
propagators, respectively. 

Once the spectral functions of the auxiliary particles 
are obtained, the physical spectral function, Pd{^) — 
— 1/nXGdiu), follows from direct convolution of the aux- 
iliary ones and the vertex functions, where G £ j(w) repre- 
sents the retarded Green function of the QD. 

It should be pointed out that the NCA equations are 
recovered from OCA when the vertex corrections are ne- 
glected in the self-energies and physical Green's functions 
bubbles. 

As already stated our aim in this work is to analyze 
the role of the crossing diagrams on equilibrium transport 
properties. In particular, on the conductance through a 
QD within the linear-response regime which is given byJ^ 



du [-/'M] A T (uj)p d {u}). (2) 



In eq. © A t (oj) = A l (lu)A r (uj)/A(uj), where A(w) = 
Al(lu) + Ar(ili) represents the total coupling between 
leads and QD, and f(ui) = 1/(1 + e^") is the Fermi func- 
tion. 



III. NUMERICAL RESULTS 

In this section, we present the results obtained for spec- 
tral density and conductance as a function of tempera- 



ture and position of the QD level within, both, NCA 
and OCA schemes. When solving the self-consistency 
equations, we consider the case of symmetric couplings, 
Al = An, and the hybridization functions to be step 
functions A u (u) = A U 6(D — \u\) with bandwidth D sev- 
eral times larger than A u . From now on, we choose the 
total hybridization A = 1 as our unit of energy. 

The numerical procedure used to solve the NCA and 
OCA equations, follows the method proposed by Kroha 
et al*£- to describe the spectral densities of the auxiliary 
and physical particles. The procedure guarantees the res- 
olution of the sets of integral equations for the auxiliary 
particles' self-energies to a high degree of accuracy down 
to very low temperatures. 

As a first step, we obtain the physical spectral function 
for the symmetric Anderson model (Ed = —U/2) within 
NCA and OCA. Fig. H shows the results for E d = -6 
at the same temperature T = O.OIT^- in units of the 
corresponding T K , being Tg CA = 0.00012 and T° CA = 
0.007. For the definition of Kondo temperature, Tk, we 
adopt the temperature for which G(Tk) = Go/2 with 
Go = 2^-. 
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FIG. 2. Physical spectral function of the QD for the symmet- 
ric case, Ed = —U/2 = —6 and D = 10. The insets show a 
zoom of the region near the Fermi energy uj = 0. 

As it was shown previously by Haule et al.^ from 
Fig. [51 it is clear that the charge transfer peaks, lo- 
cated around Ed and Ed + U, are insensitive to vertex 
corrections. On the other hand, the width of the Kondo 
resonance peak (left inset of Fig. [5]), which also repre- 
sents the Kondo temperature, increases orders of magni- 
tude within OCA as compared to NCA. The larger value 
of Tk provided by OCA can be traced back to the addi- 
tional spin-flip processes incorporated by the vertex cor- 
rections. In addition, when scaled by the corresponding 
Tk, the NCA resonance is broader than the OCA one for 
energies beyond Tk (right inset of Fig. [5]) . 

We are interested here in the modifications induced on 
the conductance by the enhancement of the Kondo scale 
and the different shape of the Kondo resonance when 
scaled by the corresponding Tk- In Fig. [3] we present 
the conductance as a function of the energy level Ed, 
which is proportional to gate voltage, in NCA and OCA 
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schemes at a high temperature, T = 0.8. It is also shown 
in Fig. [3]the NCA occupation of the QD as a function of 
Ed- We observe that there is no difference between NCA 
and OCA occupations. This results from the fact that 
this magnitude is a static one, i.e. energy-integrated, 
and therefore weakly dependent on temperature. The 
value of the occupation number indicates the different 
regimes of the QD. In particular, (n) ~ 1 stands for the 
Kondo regime. As it is seen, at this temperature value, 
there is no formation of the Kondo resonance within NCA 
nor within OCA. The conductance presents two separate 
peaks at the gate voltages for which the charge transfer 
peaks traverse the Fermi level. Due to the absence of the 
Kondo peak at this temperature (see Fig. 1 of RefJ^) 
there is no qualitative difference between NCA and OCA 
conductances, further reinforcing what we pointed in the 
spectral density of Fig. [5] 




FIG. 3. Conductance in the linear-response regime and occu- 
pation (dotted curve) of QD at high temperature, T = 0.8, 
as a function of the energy position Ed- Parameters: U = 32, 
D = 43. 

When the temperature decreases, due to the formation 
of the Kondo resonance at the Fermi level, the conduc- 
tance valley showed in Fig. [3] tends to raise as soon as 
the occupation in the QD is close to one. This means 
within the Kondo regime. 

For low enough temperatures as compared to the 
Kondo one, the Kondo effect is fully developed and the 
Kondo peak of the spectral density reaches the limit im- 
posed by Friedel's sum rule (p(0) ~ 1/ttA). As a con- 
sequence, the conductance valley evolves into a plateau 
at Go- This feature of the conductance is shown in the 
left panel of Fig. [4] for T = 0.005 within OCA and is in 
agreement with the behavior found by solving the spec- 
tral density of the QD using the slave bosons approach 
and NRC±£ 

However, as it is clear from Fig. [2j NCA results differ 
from OCA ones at low temperatures. Indeed, as we show 
in the left panel of Fig. [H due to the absence of vertex 
corrections the NCA conductance still presents a valley 
or dip for the same temperature as compared to OCA. 
To obtain a similar plateau within NCA, using the same 



set of parameters, it would be necessary to go to temper- 
atures several orders of magnitude smaller than within 
OCA. When compared at the same temperature in units 
of each Tk for the symmetric case, right panel of Fig. 
21 we observe a similar behavior of both schemes near 
this Ed value. In this case, we present the calculations at 
T = 0.06 Tk, for which the unitary limit is not exceeded 
within the Kondo regime (symmetric case). However, 
when some charge fluctuations are allowed, we obtain a 
larger violation of the unitary limit for the NCA con- 
ductance as compared to the OCA one. This represents 
an important improvement of the OCA performance over 
the NCA one. 
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FIG. 4. Conductance through the QD at low temperatures 
as a function of the energy position Ed- Parameters: U — 8, 
D = 10. 

As we mentioned previously, when comparing the spec- 
tral function in units of Tk (right inset of Fig. [5]) and 
the conductance as a function of Ed for the same tem- 
perature in units of Tk (right panel of Fig. H]), NCA 
and OCA exhibit different scaling features. Thereafter, 
in what follows, we focus on the study of the universal 
scaling properties of the conductance as a function of 
temperature. NCA results for the universal behavior of 
the conductance as a function of T/Tk, were previously 
analyzed for different values of the energy level Ed ^ Also 
for different values of Coulomb repulsion U close to the 
limit U — >• ooi£. In Fig. [5] we present the results for the 
conductance as a function of T /Tk for the symmetric 
case within the Kondo regime and at very low tempera- 
tures. In addition, we include the empirical formula for 
G(T) extracted from calculations of the numerical renor- 
malization group, NRG, for spin 1/2 



G E (T) 



G(0) 



[1 + (2V« - \){T/T K f 



(3) 



with s = 0.22^ 

Due to the overestimation of Friedel's sum rule within 
both diagrammatic schemes (see Fig. [2]) in Fig. [5] 
we normalize the conductance to G(0) = G(To), where 
To ~ O.OITk is the lowest temperature that can be 
reached without obtaining unphysical results. Below this 
temperature, we observe that the conductance still in- 
creases when a saturation is expected and therefore it 
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FIG. 5. Conductance as a function of T/Tk for the symmetric 
Anderson model. Parameters: Ed — 6, U — 12, D = 10. The 
solid, dashed, and dot-dashed lines stand for NRG, OCA, 
and NCA techniques, respectively. The inset shows the OCA 
conductance for the symmetric case, U = 12, and far from 
this, U = 18 (dot-dot-dashed line). 



can be considered as the limit of accuracy for calculating 
the conductance. 

We find that the NCA and OCA conductances follow 
different scaling behaviors in the whole range of tem- 
peratures. The OCA conductance is closer to the ex- 
act prediction of the NRG than NCA for both, low and 
high temperatures. In particular, an important improve- 
ment is obtained for temperatures larger than Tr. It 
must be noticed that vertex corrections not only change 
the Kondo scale but also the whole temperature depen- 
dence. However, it can be observed from Fig. [5] that the 
vertex functions introduced by OCA are not sufficient to 
recover the NRG prediction. As a last remark, it must 
be pointed out that while the NCA in its U — > oo limit 
follows the empirical law given by eq. we observe 

that the finite U versions presented here fail for the sym- 
metric case. At this point, a set of self-consistent vertex 
equations is still needed. The inset of Fig. [5] shows the 
conductance for the symmetric case with U = 12 and far 
from this, U — 18. We notice that for larger values of 
U, the extra diagrams not included within OCA give a 



small contribution and as a consequence the conductance 
agrees very well with the NRG prediction. 

IV. CONCLUSIONS 

In this work we obtain the conductance as a func- 
tion of both, temperature and gate voltage, for a quan- 
tum dot modeled by the Anderson impurity Hamiltonian. 
We calculate the spectral density of the quantum dot 
by using the finite Coulomb repulsion Non Crossing and 
One Crossing Approximations. The comparison of both 
schemes let us conclude about the role of vertex correc- 
tions when calculating transport properties. 

At high temperatures, as compared with the hybridiza- 
tion strength, there is no qualitative difference between 
NCA and OCA conductances as a function of gate volt- 
age as well as between the occupations of the dot. On 
the contrary, at low temperatures, the OCA conductance 
as a function of gate voltage displays a plateau, while the 
NCA one still shows a dip in the Kondo regime. This re- 
sults from the underestimated Kondo scale within NCA. 

When studying the scaling properties, we find that 
the additional processes incorporated by the vertex func- 
tions modify the functional temperature dependence of 
the conductance in addition to the Kondo temperature. 
Thus, NCA and OCA show different scaling behaviors in 
the whole range of temperatures. 

Finally, when compared with the NRG results, the 
OCA conductance as a function of temperature is more 
reliable than the NCA one. For large values of the 
Coulomb repulsion, away from the symmetric case, the 
OCA and NRG conductances agree very well. However, 
close to the symmetric case, in order to recover the NRG 
prediction, it is still necessary to go beyond OCA correc- 
tions. 
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